Multistep kinetic self-assembly of DNA-coated colloids 
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Self-assembly is traditionally described as the 
process through which an initially disordered 
system relaxes towards an equilibrium ordered 
phase only driven by local interactions between 
its building blocks. However, This definition is 
too restrictive. Nature itself provides examples of 
amorphous, yet functional, materials assembled 
upon kinetically arresting the pathway towards 
the ground state [T|-[5]. 

Kinetic self-assembly is intrinsically more flexible 
and reliable than its equilibrium counterpart, 
allowing control over the morphology of the final 
phase by tuning both the interactions and the 
thermodynamic pathway leading to kinetic arrest 

m- 

Here we propose strategies to direct the gela- 
tion of two-component colloidal mixtures by 
sequentially activating selective interspecies 
and intra-species interactions. We investigate 
morphological changes in the structure of the ar- 
rested phases by means of event driven molecular 
dynamics (MD) simulations and experimentally 
using DNA-coated colloids (DNACCs). Our 
approach can be exploited to finely tune the 
morphology of multicomponent nano- or micro- 
porous materials with possible applications in 
hybrid photovoltaics, photonics and drug deliv- 
ery. 

Thanks to the specificity of Watson-Crick base paring, 
DNACCs are considered as an ideal test bench for 
equiUbrium self-assembly. Due to the strength and 
the sharp thermal activation of the DNA-mediated 
interactions [313] j DNACCs reveal a strong tendency to 
kinetically arrest into amorphous phases, precursors of 
crystalline ground-states. Crystallization is only reliably 
achieved using nanoscale particles, p^l[T6] for which 
interactions are weaker and their thermal activation 
smoother, whereas scaling up to slightly larger building 
blocks turns out to be difficult [T71 [TB] . 
The properties that make DNACCs unsuitable for 
crystallization experiments make them ideal for the 
purpose of controlling their kinetic arrest in amorphous 
structures. We exploit both the selectivity and the sharp 
thermal activation of hybridization interactions to design 



binary mixtures of colloids in which interactions are 
activated at different stages of the quenching procedure. 
One of the species aggregates before the other, providing 
a template or a scaffold for the aggregation of the 
second. 

In contrast to equilibrium self-assembly of DNACCs, 
only possible for nanoscale objects, our approach can be 
applied to the whole range of colloidal sizes. 
In this letter we discuss two schemes that exemplify 
the potentiality of our strategy. Similar criteria can be 
applied to design arbitrarily complex systems with more 
than two components, multiple aggregation stages or 
modified interaction schemes. 

In the final section, we describe some of the expected 
technological applications. 

The first binary system is composed of species a 
and 13. The pathway to kinetic arrest is summarized in 
Figure [l] a as obtained with MD simulations. Initially 
the two-component mixture is equilibrated at high 
temperature, where it behaves as a one-component 
gas of hard spheres. A quench is then performed by 
suddenly lowering the temperature. At this stage, 
a short-range square-well attraction dominates be- 
tween a colloids, which form a gel coexisting with 
the P gas. Finally, an identical square-well attraction 
is activated between /3 colloids, which gel in the con- 
fined environment imposed by the preexisting a network. 

The two-step kinetic arrest designed in silico is tested 
experimentally using a binary mixture of fluorescently 
labeled DNA-coated microspheres with diameter of 
500 nm. 

As sketched in Figure [TJb, colloids a are functionalised 
with a symmetric mixture (ratio 1:1) of double- 
stranded DNA (dsDNA) sequences terminated with 
single-stranded DNA (ssDNA) sticky ends A and A', 
complementary to each other. The use of dsDNA rigid 
spacers allows the sticky ends to explore large volumes 
around the grafting site and limits the entropic cost 
for hybridization [7H3 [HI [SHll]- Details about the 
preparation of DNA coatings are provided in the Sup- 
plementary Information (SI). Complementarity between 
ssDNA grafted on the same particle results in the 
formation of intra-coUoid loops and inter-colloid bridges 
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FIG. 1. Self-assembly pathway and experimental implementation of the two binary mixtures, a, From left to 
right, snapshots from MD simulations of the two-step gelation process of a — /3 mixture. Initially, colloidal interactions are 
repulsive and the system equilibrates into a uniform gas phase. The a — a attraction is switched on, and gelation occurs. Finally 
the P — P attraction is activated and the P phase aggregates in the confined environment imposed by the a template-gel. b, 
Experimental implementation of the a and j3 particles by means of DNACCs. a — a and P — P hybridization interactions are 
activated at temperatures Tc, > Tp such that the two-step aggregation pathway can be induced upon quenching from Tin > Ta 
to Ta > Tmid > TJg and finally to Tfln < T/3. c. From left to right, snapshots from MD simulations of a core-shell gelation of 
7 — 5 mixture. Initially the 7 — 7 attraction is switched on, and 7 particles form a gel from the initially uniform gas phase. 
Upon activation of the 'y — 5 attraction, S particles coat the 7 gel as a single layer. 7 — 7 interactions are repulsive at every 
stage, d, Experimental implementation of the y — 5 mixture with DNACCs. Colloids 7 are functionalised with B and B' sticky 
ends, which produce interspecies attraction for T < Tj. Colloids 5 are coated only with strand Bs, complementary to part of 
B' (marked in red), such that a 7 — 5 attraction arises for T <Ts < T^. The DNA sequences are grafted to the surface of the 
colloids via biotin-streptavidin bonds. To reduce the aggregation temperature of species a and 7 a fraction of the grafting sites 
are occupied by inert polyethylene-glycol (PEG) polymers. 



[T71 [T^MI]- The latter are responsible for interspecies 
attraction between a colloids, activated by quenching 
below the aggregation temperature Ta ~ 43° C. 
Similarly to a, colloids /3 are coated with a symmetric 
mixture of complementary Bg and Bg strands. The 
aggregation between /? particles is however triggered at 
temperature Tp « 33°C< T„. 

The two-step gelation is replicated by quenching the 
a — /3 mixture from the initial gas phase, allowed to equi- 
librate at Tin = 55° C, to an intermediate temperature 
Tmid = 40°C, with < Tmid < T^. At this stage, a 
colloids undergo gelation, whereas /3 colloids remain in a 
gas phase. After a percolating a network is formed, we 
further decrease the temperature to Tfjn < T/j to induce 
gelation of the /3 species. Epifluorescence microscopy 
snapshots of the three stages are shown in Figure [2] a. 
We construct samples in a slab geometry to ensure high 



quality imaging, details are provided in the Methods 
section and SI. 

We monitor the evolution of both the species by mea- 
suring the quantity M(t) = max [Sq{t)/Sq{t = 0)], where 
Sq{t) — (/q/_q)/ J I{x,y)dxdy. The intensity /q is the 
Fourier transform of an epifluorescence snapshot /(x, y) 
and (...) indicates a radial average. The value Sq{t = 0) 
is measured for the gas phase at T = Tin. The evolution 
of M{t) for the a — (3 mixture, shown in Figure [2] b, de- 
picts the two-step gelation. 

We obtain morphological information about the arrested 
structures from confocal imaging. We produce 3D re- 
constructions by processing confocal stacks with Gaus- 
sian filtering and thresholding algorithms (SI). Quantita- 
tive characterisation is carried out by means of a chord- 
distribution analysis [53J[21]. As shown in Figure S4, the 
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FIG. 2. Two-step gelation with DNA-coated colloids, a, Epifluorescence snapshots of a (top, red frame) and /3 (bottom, 
green frame) colloids at different stages of the quenching procedure for a sample with (j>ci = (j^P ~ 5%. The scale bar is equal 
to 20 /^m. b, Evolution of temperature (right axis) and M{t) — max[Sq{t) / Sq{t = 0)] for the two components (left axis) for a 
sample with 0^ ~ 5% and <j!>/3 ~ 10%. c, Confocal fluorescence images of three arrested samples and corresponding top-view of 
3D volume reconstructions with ~ 5% and, from left to right, (^fj ~ 5, 7.5, 10%. The scale bar is equal to 50 /^m. d. Chord 
distributions for the j3 phase of samples shown in panel c (o) and exponential fits (solid lines). In the inset: persistence length 
A evaluated from chord distributions of a and /3 components for samples shown in panel c. 



chords are measured as the lengths of the intersections 
between straight hnes drawn on the thresholded images 
and regions belonging to the aggregates. Chord length 
histograms /(r) exhibit the exponential decay typical of 
patterns originating from arrested phase separation. The 
decay length A is defined as the persistence length of the 
structures [ST and provides a measure of the length-scale 
of the kinetic arrest |23[ . 

In Figure [2] c, we show coexisting a — /3 gels with a con- 
stant volume fraction for the a phase (f)a ~ 5% and in- 
creasing « 5, 7.5, 10%. The morphology of the /3 
phase is influenced by the preexisting a gel. This is ev- 
ident in the case of symmetric composition = 0^, 
for which, due to the confinement imposed by the a gel, 
the (3 species aggregates into disconnected clusters. Such 
morphological asymmetry is confirmed by the chord- 
distribution length-scales Aq, > shown in the inset 
of Figure [2] d and is a direct consequence of the two-step 
kinetic arrest. In fact, a — /3 mixtures with = Tp 
form coexisting gels with identical morphology |22) (see 
SI and Figure S5). Xp grows upon increasing (j)/^ to 7.5 
and 10%. For the sample with = 2(^c< we observe a 
significant drop in Aq., indicating that a high concentra- 
tion of /3 colloids in the gas phase inhibits the gelation of 



the a component. 

Besides the confining effect of the a gel, the morphology 
of the /3 phase can be influenced by nonspecific a — /3 at- 
tractions that might result in the trapping of /3 particles 
within the a aggregates. By introducing controlled in- 
terspecies interactions mediated by DNA, we could take 
advantage of this effect and further manipulate the struc- 
tures. 

This kinetic route can be exploited to produce struc- 
tures characterised by an even higher level of complexity. 
We exemplify this by presenting a second binary mixture 
composed of colloids 7 and 6. The two-step arresting 
sequence is illustrated in Figure [l] c as obtained from 
MD simulations. After equilibrating the system at 
high temperature, we trigger the gelation of the 7 
component by activating the 7 — 7 attraction, analogous 
to the first gelation step of the a — /3 mixture. For the 
second step, we activate an interspecies 7 — ^ attraction, 
maintaining the S — S interaction repulsive. This results 
in the formation of a single layer of S particles coating 
the percolating 7 gel, and producing a core-shell gel 
structure. 
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FIG. 3. Core-shell gels with DNA-coated colloids, a, Epifluorescence snapshots of 7 (top, red frame) and /3 (bottom, 
blue frame) colloids at different stages of the quenching procedure for a sample with (jt-y = it's ~ 5%. The scale bar is equal to 
10 /im. In Supplementary Video 2 we show the evolution of the two-step gelation for a sample with 0^, « 5% and (pg ~ 7.5%. 
b. Evolution of temperature (right axis) and of M{t) = max [Sq{t) / Sq{t = 0)] for the two components (left axis) for a sample 
with (fi-f — cj>s ~ 5%. c, Confocal image and corresponding 3D reconstruction for a sample with (j)~t ~ 4>& ~ 5%. The scale bar is 
equal to 50/im. d. Chord distribution calculated for the sample in panel c. Blue o indicate the 5 phase, red o the 7 phase and 
grey o the total network obtained by adding up the & and 7 images. The lines are exponential fits, e. Higher magnification 
confocal images showing samples with constant (f)-^ = 5% and, from left to right, 05 = 5, 7.5, 10%. For the last sample we 
show a time-averaged image to distinguish the 5 gas phase from the 5 shell surrounding the 7 core gel. f. Fraction /coat of 5 
colloids adhering to the 7 gel as a function of calculated with MD simulations. 



Similarly to a, 7 DNACCs are functionalised with 
a symmetric mixture of complementary sticky ends 
B and B' , which cause aggregation at temperature 
T < ^ 43° C. Colloids 5 are functionalised only 
with Bg sticky ends, complementary to a subsection 
of B' (Figure [T] d). Upon quenching to temperatures 
below Ts « 33° C a specific •y — 6 attraction is activated, 
which triggers the formation of the coating. Since Bg is 
not self-complementary, the S — S interaction remains 
repulsive regardless of temperature. 

The two-step formation of core-shell gels is demon- 
strated by epifluorescence images in Figure [3] a and 
by M{t) in Figure |3] b, both relative to a sample with 
4>-f = <t>5 ~ 5%. In Figure [3] c we show a confocal image 
with its 3D reconstruction. A and a zoomed confocal 
stack is shown in Supplementary Video 1. Also for the 
core-shell structures, the chord distribution provides 
quantitative morphological information. We find a 
persistence length of the coating phase Xs ~ 570 nm, as 
expected, similar to the diameter of a single colloid. The 
persistence length of the total gels \^+s is consistently 
larger than A^. 

To further demonstrate that the 5 — 5 attraction is 



negligible and 5 colloids only form a single layer coating 
the surface of 7 gel, we produce samples with increasing 
and constant (p-^ 5%. High magnification confocal 
images are shown in Figure |3]e, full sized ones in Figure 
S6. For the case of 4)s = 5%, the amount of 5 colloids in 
the gas phase is similar to the amount of 7 colloids. Both 
are compatible with a gas-gel coexistence. Increasing <j)^ 
to K, 7.5%, we observe a clear growth of the number of 
free 5 colloids, indicating that the surface of the 7 gel 
is saturated and further adhesion of the 5 particles is 
hindered. This effect is more evident for samples with 
(ps — 10%, for which a time-averaged image is necessary 
to distinguish between the immobilized 5 particles 
making up the shell and the abundant 5 gas. 
We quantify the saturation effect with MD simulations 
by measuring the fraction /coat of 5 colloids adhering to 
the 7 gel as a function of (jis/cji-yj as shown in Figure [3] f. 
In agreement with the experiments, the saturation is 
observed to occur for 1 < (j)s/ (j)-^ < 1.25. 

Multi-step kinetic self-assembly enables further mor- 
phological manipulation by changing the timescale of 
the quench process. We can process the same a — 13 and 
7 — (5 mixtures with a very fast, one-step, quench instead 
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of the two-step one such that the second component is 
activated before the gelation of the first. For the a — /3 
mixture, this removes the morphological asymmetry, 
while the 7 — 5 mixture aggregates in a single binary 
gel. All the morphologies between this limiting situation 
and the ones obtained the with two-step quench are in 
principle accessible and will originate materials with 
different tunable properties. 

Colloidal structures assembled via DNA can be made 
permanent |20j . enabling a direct application of our 
strategy to the self-assembly of functional materials. 
For example, spongy phases generated in nature by ar- 
rested phase separation, can produce structural coloring 
[2H5]. Similar materials can be engineered based on 
our approach, with the additional possibility of tuning 
refractive index contrast by having multiple components. 
In random lasers, amorphous assemblies of colloids are 
immersed in a medium that provides optical gain to 
ensure the multiple scattering required to reach the 
lasing threshold [53] . Resonant emission can be obtained 
by using randomly packed monodisperse spheres as 
scatterers [Mj . Similarly, we suggest the use of colloidal 
gels with tunable photonic response. Multiple compo- 
nents would also allow decoupling between the spatial 
distribution of the dyes providing gain, which can be 
incorporated in latex colloids [17], and the scatterers. 
In hybrid photovoltaics, controlled phase separation 
is required between the conjugate polymer acting as 
light absorber/donor and inorganic particles acting as 
acceptors |28|. Schemes analogous to our core-shell 
self-assembly could be exploited to build inorganic 
semiconductor scaffolds of different morphology and 
coat them with organic particles. 

Furthermore, the reversibility of the DNA interactions 
can be exploited to engineer responsive materials for 
sensing and drug release. Let us consider the core-shell 
gels. The scaffold phase can be made permanent and 
used as a means of containing and delivering controlled 
quantities of the coating colloidal species. The latter 
can then be released by an external stimulus, such as 
change in temperature or pH. 

We conclude by stressing that our strategy is not limited 
to the use of DNA. Supramolecular chemistry provides 
a wide range of tags with selective interactions [25], the 
functionality of which is not limited to physiological 
environments as for DNA. 

Summarizing, we introduce an alternative means of 
controlling the self-assembly of colloidal particles into 
amorphous materials. Contrary to equilibrium self- 
assembly strategies, in which the target configuration 
only depends on the interactions between the building 
blocks, our approach relies on both the selectivity and 
the temperature dependence of the base-paring interac- 
tions and exploits multi-step kinetic arrest to tune the 



morphology of the final configuration. 
We describe two model systems to demonstrate the po- 
tential and flexibility of our scheme. The first is a binary 
mixture undergoing a two-step gelation, with the second 
aggregating phase forced to gel in the confined environ- 
ment produced by the first aggregating phase. The sec- 
ond system self-assembles in core-shell colloidal gels, with 
the second component coating the network firstly formed, 
which acts as a scaffold. 

Multi-step kinetic self-assembly of mesoscopic building 
blocks mediated by DNA specific interactions holds great 
promise as general approach for reliable self-assembly of 
mesoporous materials. 

Methods 

Experimental methods 

DNA (Integrated DNA technology) is grafted to the surface of 
streptavidin-coated colloids via biotin-streptavidin linkage. The 
DNA constructs are structured as follows: Biotin-5'-TTTTT- 
dsDNA spacer-TTTTT-ssDNA sticky end-3'. The dsDNA spacer 
has a length of 60 base pairs. The 5-thymine sequences are added 
to confer more flexibility to the constructs. The melting tempera- 
tures of the free sticky ends in solution are: B — B': 37.5°C, A — A': 
37.8°C, Bs - B'g{A'): 30.2°C, calculated at 50mM ionic strength. 
We use two batches of green and red fluorescent polystyrene col- 
loids (Microparticles GmbH, Berlin). The experiments are carried 
out in 10-mM Tris-HCl pH 8 + 1-mM EDTA bufl'er with the addi- 
tion of 50-mM NaCl. The density of the solution is raised to 1.05 
g/cm'^ to match that of polystyrene. Due to the high refractive 
index of polystyrene and the incompatibility of DNA with organic 
solvents, refractive index matching limitations exist. To allow good 
quality optical imaging of opaque samples, glass sample chambers 
are built in a slab geometry, with width of 1.8 cm and thickness of 
15 ^m. The chambers are permanently sealed to avoid evaporation. 
Further details are provided in SI. Heating/quenching procedures 
are carried out while imaging with a Nikon Eclipse Ti-E inverted 
microscope using Plan Apo VC 20x/0.75 or Plan Fluor 40x/0.75 
dry objectives. Confocal imaging is performed on a Leica TCS 
SP5 microscope equipped with an HCX PL APO CS 100 X 1.4 oil 
immersion objective. 

Computer simulations 

Event-driven molecular dynamics simulations are performed in 
cubic boxes with periodic boundary conditions. The colloids are 
modeled by means of a square-well potential, with a hard-sphere 
core and with a well whose depth is either set equal to eij = —eg 
when the i and j species are attractive, or to eij = when they 
are considered as repulsive. The simulations involve up to 2 X 10** 
particles. The colloids have equal diameter cr = 1 and equal mass 
m = 1. The potential is chosen to be short-ranged, i.e. the 3% of 
the diameter. We used reduced units, kg = 1 and eo = 1. The 
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quenches are performed by lowering the temperature from T = 100 
to the final value T = 0.05, well below the known critical temper- 
ature Tc = 0.3 for a single-component system |30) . As the system 
falls out of equilibrium, averages over 10 independent realizations 
of the initial conditions are considered instead of time averages. 
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